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Effects of the site dilution on spin-gapped Heisenberg antiferromagnets with 5 = 1/2 and S = 1 
on a square lattice are investigated by means of the quantum Monte Carlo method. It is found that 
effective magnetic moments induced around the diluted sites exhibit the antiferromagnetic long- 
range order in the medium of spin-singlet pairs. Their microscopic structure is examined in detail 
and important roles of the higher dimensionality than one on the phenomenon are discussed. 
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The nonmagnetic-impurity-induced, or more specifi- 
cally, the site-dilution-induced antiferromagnetic long- 
range order (AF-LRO) in quasi-one-dimensional (Q1D) 
quantum antiferromagnets has attracted many re- 
searchers more than a decade jjj. Typically it has been 
observed in the first inorganic spin-Peierls (SP) com- 
pound CuGeC>3 Q by the substitution of nonmagnetic 
Zn for magnetic Cu with 5 — 1/2 spin. In the induced 
AF state with sufficiently weak dilution the lattice dimer- 
ization is preserved ||. This strongly suggests that at 
least this site-dilution-induced transition is dominantly 
of a magnetic origin, though for thorough understand- 
ing of the SP transition in the system other ingredi- 
ents such as the magneto-elastic effect Qj is indispens- 
able. Another example of the site-dilution-induced AF- 
LRO has been observed in a Q1D Haldane compound 
Pb(Nii_ a ,Mg 2 ,) 2 V 2 8 where Ni carries 5 = 1 spin |g. 

A most plausible scenario for the site-dilution-induced 
AF-LRO so far discussed is as follows ||||]. The 
mother compound (x — 0) has a spin-gapped ground 
state such as the dimer singlet state in CuGeOa and the 
Haldane state in PbNi2V2 0s. When x > 0, spins on the 
sites neighboring the diluted sites organize a spin clus- 
ter which is a close resemblance of the edge state in the 
Haldane chain E3], Here we call it as a whole an 'ef- 
fective spirt. Between the effective spins there exist the 
effective interactions mediated by spins which are in the 
singlet state. The averaged magnitude of the effective 
interactions is exponentially small as a function of the 
concentration of nonmagnetic impurities. Although they 
are either ferromagnetic or AF depending on relative po- 
sitions of the effective spins, the staggered nature with 
respect to the original lattice is completely preserved, and 
so the effective spins are antiferromagnetically ordered at 
T = in more than two dimensions. The very last part 
of the above scenario, i.e., the effect of the higher dimen- 
sionality, has been mostly treated by the mean-filed-like 
arguments except for a few works . 

In the present work we have extensively performed the 
quantum Monte Carlo (QMC) analysis on the site-diluted 
2D Heisenberg antiferromagnets on a square lattice, and 
have reached the same scenario, but with evidences which 



support it more directly than the previous works. One 
of them is the microscopic spatial information on the ef- 
fective spins induced by the dilution. We also present 
the result consistent with the scenario in the 2D coupled 
Haldane chain system with 5 = 1. 

The quantum spin system of our present interest is 
described by the Hamiltonian 

1~i = ^2i,j^2i+l,j^2i,j ' S2i+lj 

+ a £2i+l ,j^2i+2,j^2i+l.j ' ^2i+2.j (1) 

ij 

i,j 

where 1 and a (0 < a < 1) are the intrachain alternating 
coupling constants, J' (> 0) is the interchain AF coupling 
constant and Sjj is the quantum spin operator at site 
(i, j). Randomly quenched magnetic occupation factors 
independently take either 1 or with probability 
1 — x and x, respectively, where x is the concentration of 
the diluted sites. 

In the present work the QMC simulations with the 
continuous-imaginary-time loop algorithm p"l"| ] is carried 
out on the 5=1/2 and 5 = 1 systems on L x L square 
lattices with the periodic boundary condition. For di- 
luted systems with x > the number of samples aver- 
aged over is 100 ~ 700 depending on x. For each samples, 
10 4 Monte Carlo steps (MCS) are spent for measurement 
after 10 3 MCS for thermalization. The staggered mag- 
netization, M s (x), at zero temperature is evaluated by 

M s 2 (x) = lim lim 3S *( L > T > X ) ( 2 ) 

where 

4(1,1,^13-1)1^1^) (3) 

id 

is the static staggered structure factor. The bracket (■ • •) 
in Eq.([|) denotes both the thermal and random averages. 
The value of 5 S (L, T, x) converges to its zero-temperature 
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FIG. 1. Concentration dependence of the staggered mag- 
netization at zero temperature for systems with a = J' = 1, 
0.6, 0.55, 0.5, and 0.4. The a — J' phase diagram of the pure 
systems is shown in the inset. The data for a = J' = 1 are 
results in ref. [13] with :r c i(— 0.407254) being the percolation 
threshold. All the lines are guides to eyes. 



value at T lower than either a gap due to the finiteness of 
the system or the intrinsic spin gap. Thus S s (L,T,x) at 
low temperatures where its T-dependence becomes not 
discernible within the error bars is taken as an estimate 
of S s (L,0,x). 

First we determine the ground-state phase diagram of 
the pure systems on the a — J' plane. Its details will be 
discussed elsewhere |l2| . Here we simply show the result 
for the 5=1/2 system in the inset of Fig. 1, and mention 
that, except for the critical point (ex, J') = (1,0), the 
quantum dimer-AF transition we obtain is ascertained 
to belong to the same universality class as that of the 3D 
classical Heisenberg model. 

Now let us concentrate on the site-dilution-induced 
AF-LRO in the S = 1/2 systems with a = J' indi- 
cated in the inset of Fig. 1. Expecting that the fun- 
damental aspects of the phenomenon do not crucially 
depend on the values of a and J', we choose those on 
the line a — J'. It connects the system with a = J' = 0, 
which consists of a set of independent dimers, and the 
isotropic antiferromagnet with a = J' = 1. Along this 
line, the quantum dimer-AF phase transition occurs at 
a = a c = 0.52337(3) @. The dilution effects on the iso- 
toropic system was investigated in details in the previous 
works |l3|jl4f| : as shown in Fig. 1, the site-dilution sim- 
ply reduces M s which exists already in the pure system. 
For the systems with a = J' = 0.5 and 0.4, on the other 
hand, M s , which is zero for x — 0, becomes finite even 
with 1 or 2 % of site dilution. This is nothing but the 
site-dilution-induced AF-LRO of our present interest. 

The values of M s in Fig. 1 are extracted from the 
L-dependence of S s (L,0,x)/L 2 as shown in Fig. 2. 
The data for x = are well fitted to S s (L,0,x)/L 2 cx 
(1 — exp(— L/£ p ))/L 2 , where £ p is the correlation length. 
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FIG. 2. System-size dependence of S S (L, 0, x)j L 2 in the 
case of a = J' = 0.5. Dashed lines are obtained by the 
least-squares fitting for the largest three system sizes for each 
x. The extrapolated values are denoted by solid symbols. 



This form is derived by the modified spin-wave theory 
for a spin-gapped state. For diluted systems with x > 0, 
on the other hand, M s is evaluated by fitting the data to 
S s (L,0,x)/L 2 ~ M s 2 /3 + a/L which is obtained, for ex- 
ample, by the linear spin-wave theory on the state with 
a finite M s . 

The above-mentioned correlation length of the pure 
system, £ p , is given by £ p = (C p ^) 1/2 where and 
are the anisotropic correlation lengths. They are esti- 
mated from the dynamic staggered correlation function 
by using the second moment method. Similarly the spin 
gap, A p , is evaluated from behavior of the function on 
the imaginary-time axis. We obtain |lj] £ p = 3.0089(9), 
^ = 2.2097(6) and A p = 0.32255(3) for a = J' = 0.4 
and = 11.998(9), f» = 9.312(10) and A p = 0.0918(1) 
for a = J' = 0.5. The nearer to the transition point a 
system is, the smaller is A p and the larger is £ p . 

In order to gain insights into nature of the site-dilution- 
induced AF-LRO obtained above, we investigate the local 
static staggered structure factor defined by 



5(i) = £(-i)l 



\(SfS!) 



(4) 



Figure 3 shows the real-space distribution of S(i) for 
a = J' = 0.5 on a 64 x 64 lattice from which 30 spins are 
randomly removed. The temperature, T — 0.001, can be 
regarded as zero temperature for the system with these 
parameters. The points protruding downwards to zero 
are on the diluted sites. The peaks of S(i) appear on the 
sites connected by a strong bond 1 with these diluted 
sites. Naturally they are considered to be the contribu- 
tion of the spins which get rid of the singlet pairing by the 
dilution. The peaks exhibit a significant spatial extent 
whose size is given by £ p and £j( evaluated above. 

From the analogy of the ID spin-gapped state g] , the 
exchange coupling between the two effective spins cen- 
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FIG. 3. 3D plot of the local static staggered structure fac- 
tor in the fixed configuration of 30 impurities on a 64x64 
lattice for a = J' = 0.5 at T = 0.001. Contours are shown in 
the bottom. 
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FIG. 5. Temperature dependence of the effective Curie con- 
stant for a — J' — 0.4. The broken lines show \T = x/A for 
x = 0.04, 0.02, and 0.01 from top. 
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FIG. 4. 3D plot of the staggered correlation function from 
the peak of an induced magnetic moment on site (=(31,31)) 
for the same dilution configuration as Fig. 3. 



tered at sites m and n is expected to be given by J mn oc 
(-l) |rm - r " +1| exp[-/(^)- 1 / 2 ], where I = \r m - r n \ is 
the distance between the effective spins. For a suffi- 
ciently weak dilution, we can naturally expect that the 
effective Hamiltonian, H c e = Z)< mn > JmnSm • S„ well 
describes magnetic behavior associated with the site- 
dilution- induced AF-LRO. As we have already pointed 
out, J m n is considered to completely preserve the stag- 
gered nature with respect to the original lattice since 
the original interactions are all AF with no frustration 
at all. This is in fact confirmed by Fig. 4 where we 
draw the static staggered correlation function C(i,0) = 
(-l)\n-ro\(s?S§) with site (= (31,31)) being the peak 
position of an effective spin. The system is the same as 
that in Fig. 3. At any site i, C(i,0) is positive except 
for the diluted sites where it is zero. It peaks at sites 
of other effective spins with heights almost independent 
of the distance from site 0, while it stays at a relatively 
lower value on the sites far from the diluted sites where 
spins form the singlet state. Thus Fig. 4 combined with 



Fig. 3 clearly demonstrates that the site-dilution-induced 
AF-LRO is mostly carried by the effective spins. Each 
of them contributes to M s by the amount proportional 
to £p. This explains a sharper rise near x ~ of M s of 
a = J' = 0.5 than that of a = J' = 0.4 in Fig. 1. It con- 
tributes to the uniform susceptibility, on the other hand, 
by the amount |S m | = 1/2 |6| as discussed below. 

The averaged distance (I) between diluted sites is pro- 
portional to 1/ y/x. The typical coupling which scales 
the excitation energy of the spin-wave-likc modes is 

then given by (|J mn |) oc exp(— l/^/a^pCp)- It becomes 
much smaller than A p for a sufficiently small value of 
x. Thus we can naturally expect the existence of low- 
lying excitation modes well separated from the origi- 
nal triplet excited state with A p as has been in fact 
observed by the neutron inelastic experiment in the 
doped CuGe03 [| 15 R]] . Another peculiar feature of the 
present site-dilution-induced AF-LRO is that it is inho- 
mogeneous in a shorter length scale than (I). This may 
explain the experimental fact |3j that the spectrum of the 
spin- wave-like mode becomes broader for the larger wave 
number. 

By the present simulation the existence of low-lying ex- 
citation modes mentioned above is ascertained by the In- 
dependence of the effective Curie constant \F || which 
is shown in Fig. 5 for a = J 1 = 0.4. For x = 0.01 
the plateau of a height nearly equal to \F — x/A in 
unit of g 2 /i|/fcB is clearly seen around T ~ 0.05 which 
is smaller than A p by an order of magnitude. We can 
attribute the plateau to the Curie law of the effective 
spins with |S m | = 1/2 as mentioned above. They start 
to correlate at lower temperatures than (| J mn \). Already 
for x = 0.02, however, the plateau is obscured, indicat- 
ing that (| J mn |) becomes comparable in magnitude with 
A p . Qualitatively similar features have been observed 
experimentally jlTy . 

Essentially the same scenario as the S — 1/2 case is 
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FIG. 6. The staggered magnetization of the coupled Hal- 
dane chain systems with J' = 1 and 0.04. The data for 
a = J' = 1 are results in ref. [13] . All the lines are guides to 
eyes. 

expected to hold also for similar antiferromagnets with 
5 = 1. As their representative we have examined the cou- 
pled Haldane chain system with a = 1 and J' = 0.04 in 
Eq. ([j]) . The latter is a little smaller than the critical value 
J' c = 0.04365(1) for the Haldane- AF transition point @. 
Behavior of its M s {x) shown in Fig. 6 is quite similar to 
that in Fig. 1, and further confirms the above-mentioned 
scenario for the site-dilution- induced AF-LRO. 

Here it is worth noting that the interchain interaction 
J' introduces another energy scale than A p and (|J m n|) 
into the problem of the present interest. It is the coupling 
strength, J , of spins on the two opposite neighbors of the 
diluted site through the shortest interaction paths. It is 
ferromagnetic and of the order of J' 2 . For a sufficiently 
small x, (\J mn \) < J < A p holds. In such a coupled 
Haldane chain system the following successive process 
is expected to occur; effective spins with 5 — 1/2, in- 
stead of 5 = 1, induced at both sides of diluted sites 
are independently fluctuating at J C T C A p , the 
two spins are coupled to behave as one effective spin of 
5 = 1 at (\J mn \) « T « J, and then the S = 1 ef- 
fective spins begin to correlate antiferromagnetically at 
T <C (| JmnD- F° r the Haldane chain system in Fig. 6, 
A p ^> J ~ 0.0016 is satisfied. Unfortunately however, 
the expected two-plateau structure in \T is hard to be 
seen, since J is too small. In the 5=1/2 system of Fig. 
1, on the other hand, J ~ A p , and there appears one 
plateau at T <C J as already shown in Fig. 5. 

Another interesting problem on which J plays an es- 
sential role is a possible difference between the AF-LRO 's 
induced by site-dilution and bond-dilution. In the lat- 
ter case the coupling J between the two end spins of 
a diluted bond becomes AF. Therefore, if \ J\ (~ J' 2 ) 
is significantly larger than (|J mn |), the two spins may 
again form a singlet pair and cannot contribute to the 
AF-LRO. This implies a possible existence of a criti- 



cal concentration x c below which the AF-LRO does not 
come out. In contrast, there is no reason for a finite x c 
in the site-diluted case discussed above. Experimentally, 
however, no drastic difference has been observed between 
site-diluted and bond-disordered CuGeOa yet 16|. For 



quantitative understanding of the experimental resuts we 
have to take into account other ingredients than those 
in Eq. (|l|) such as the the next-nearest-neighboring intra- 
chain interactions |ll| . This is beyond the scope of the 
present work. 

In summary, we have argued based on the QMC simu- 
lation on the 2D AF Heisenberg models that dilution of 
spins in a system having the spin-gapped ground state 
induces effective spins which are strongly correlated. As 
a consequence there occurs the AF-LRO state which has 
low-lying excitation states well separated from the origi- 
nal triplet excited state. 

Most of numerical calculations in the present work 
have been performed on the SGI 2800 at Institute for 
Solid State Physics, University of Tokyo. The present 
work is supported by the "Research for the Future Pro- 
gram" (JSPS-RFTF97P01103) of Japan Society for the 
Promotion of Science. 
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